Electron correlations and bond-length fluctuations in copper oxide superconductors: 

electron versus hole doping 
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We investigate the nature of the electronic ground state and electron-lattice couplings for doped 
chains of CuCU plaquettes or CuC>6 octahedra. The undoped configuration implies here Cu 3d 9 
and O 2p 6 formal valence states. The results of multiconfiguration calculations on 4-plaquette (or 
4-octahedra) linear clusters indicate strong electron-lattice interactions and polaronic behavior of 
the doped particles, for both electron and hole doping. For certain phases of the oxygen-ion half- 
breathing distortions a multi-well energy landscape is predicted. Since each well is associated to 
carriers localized at different sites, the half-breathing displacements induce charge transfer along the 
chain. In the case of hole-doping, the trends found by ab initio multiconfiguration calculations on 
4-octahedra clusters are confirmed by density-matrix renormalization-group calculations for a p— d, 
extended Hubbard model with chains of few tens of CuCU plaquettes. Under the assumption of 
charge separation and the formation of 1/3-doped stripes, our results seem to support the traveling 
charge-density wave scenario proposed in some recent contributions for superconductivity in copper 
oxides. 



I. INTRODUCTION 

A remarkable property of the copper oxide compounds 
is that superconductivity can be obtained in these sys- 
tems by both hole and electron doping. Holes added to 
the insulating, antiferromagnetic (AFM) CuC>2 network 
have O 2p character, see for example [E S S S H H . 
Results of recent ab initio, multiconfiguration calcula- 
tions on clusters including several CuC>6 octahedra also 
indicate strong electron - lattice couplings [f|. It was 
found in Ref. [g that the lowest energy configurations 
imply Zhang-Rice (ZR) like, quasi-localized states 
on distorted CuC>4 plaquettes. However, O-atom half- 
breathing displacements that restore the high-symmetry 
structure with identical Cu-0 distances induce strong 
charge redistribution. The O 2p hole can be (partially) 
transferred onto a single anion to give an electronic wave- 
function with a dominant contribution from a ...-Cud 9 - 
Op 5 -Cud 9 -... configuration. The ZR like singlet can 
hop thus within the CuC>2 plane via such cP-p 5 - d 9 states 
through coupling to the oxygen vibrations. According to 
the results from [(|, the energy barrier associated with 
this double- well potential is 300-400 meV for an isolated 
O 2p hole. 

For not too high doping levels, these ZR type singlets 
are sufficiently far apart such that they do not overlap. 
It seems that at the optimal hole concentration p as 1/6 
and low temperatures, spin correlations and hole-hole 
Coulomb interactions determine within the Cu02 plane 
an unique, charge inhomogeneous arrangement of the O 
2p holes, see for example £1, HI, El El El ■ J- B. Goode- 
nough [!2 argued that at p w 1/6 the doped holes form 
alternating, hole- free and 1 /3-doped stripes of CuC>4 pla- 
quettes. Using the ab initio results from Ref. for hole- 
doped clusters, the charge and spin distribution for a 
1/3-doped stripe could be schematically represented as in 
Fig. 1(a). The squares correspond to ZR states on CuC>4 



plaquettes with shorter Cu-0 bonds along the chain |(| . 
In the undistorted structure, a quasilocalized O 2p hole 
favors ferromagnetic couplings between adjacent Cu ions, 
as initially predicted in |13L fl4L fl5l| - A three-site, d 9 -p 5 - 
d 9 spin-polaron model was used by Kochelaev et al. to 
explain the electron paramagnetic resonance (EPR) sig- 
nal at different doping levels in La2- a; Sr 3 ;Cu04 |16| . 

We speculated in that for a distribution of the oxy- 
gen holes as that shown in Fig. 1(a), the strong coupling 
between the nuclear and electronic motions, the hole-hole 
Coulomb repulsion, and the in-plane spin interactions 
would cause strong correlations among half-breathing O- 
atom vibrations on every third, ZR like Cu04 plaquette 
along the doped chains. The superconductivity might 




FIG. 1: Possible representation of the traveling charge- 
density /spin-density wave along a 1/3 hole-doped CuCU- 
Cu04~... chain, see Section III and Ref. 0. The O 2p holes 
are shown schematically with black or gray color. The ellipses 
indicate spin-singlet couplings. 
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be related to the collective, coherent motion of the O 
2p holes along the 1/3-doped chains, which is associ- 
ated with phase coherence among the kind of O-atom 
vibrations shown in Fig. l(a-c). This evokes the vibronic 
mechanism and the traveling chargc-density/spin-density 
wave (CDW/SDW) scenario suggested by Goodenough 

m 

Spatially localized bound states of two electrons were 
analyzed within a two-dimensional Holstein-Hubbard 
model on a square lattice by Proville and Aubry [l7| . 
Two electrons bound into a localized singlet state as a 
result of strong electron-lattice interactions are often re- 
ferred to as a bipolaron. It was found in Ref.0] that 
for certain values of the model parameters a "quadris- 
inglet" (QS) state similar to the ZR like state discussed 
above can be formed. The near-degeneracy between the 
QS "bipolaron" and a two-site spin-singlet bipolaron can 
lead to a sharp reduction of the bipolaron effective mass. 
Proville and Aubry argued that such near-degeneracy ef- 
fects plus quantum tunneling effects associated with a 
rather flat energy landscape should be crucial ingredients 
in the theory of "bipolaronic" superconductivity 

The fact that the ab initio calculations from Ref.0 
confirm a number of distinct effects and concepts (that 
is, ZR type states on CUO4 plaquettes pj, Kondo like in- 
teractions and ferromagnetic couplings between 3cL2 ,,2 
electrons adjacent to a quasilocalized oxygen hole (l3lll4l 
Il5j . near-degeneracy between the QS configuration and 
a two-site spin-singlet configuration and strong electron- 



lattice couplings [T3I) deduced or proposed on the basis 
of model Hamiltonians that rest on quite different ap- 
proximations and describe rather different types of inter- 
actions is very interesting. It illustrates both the power 
of and the need for first principles, explicitly correlated 
methods in the study of these materials. Such calcula- 
tions on not very large clusters are sometimes essential 
for understanding the "local" many-body physics and for 
setting up the "right" model for the extended system, as 
recently shown in the case of a mixed-valence vanadium 
oxide, NaV 2 5 Near-degeneracy and pseudo Jahn- 
Tcllcr effects involving O-hole d 1 p 5 d configurations 
similar to the d 9 -p 5 -d 9 three-site spin-polaron in the 
doped copper oxides appear to be the source of most of 
the intriguing properties of this compound. 

In the present paper we report results of ab initio 
wave-function based, embedded cluster calculations for 
the electron-doped case. We investigate the ground-state 
electron distribution and electron-lattice couplings. For 
an 1-doped-electron, 4-plaquette linear cluster, our data 
show that the lowest energy state corresponds to a con- 
figuration where the extra electron is localized on a CUO4 
plaquette with elongated Cu-0 bonds. As in the hole- 
doped system, the O-atom half-breathing displacements 
interact with the electronic degrees of freedom and may 
induce charge transfer along the CuO chain. We shall de- 
scribe various aspects of this kind of interactions and dis- 
cuss the similarities and differences with the hole-doped 
situation. 




FIG. 2: Electron-doped, 4-plaquette cluster and the phases of 
the O-ion half-breathing distortions investigated by CAS SCF 
calculations. A Cu 3d^ 2 _ y 2 ion is shown with a filled, dark- 
gray circle. From left to right, the transition metal sites are 
labeled as Cun, Cu;, Cu r , and Cu rr . Ligands that are dis- 
placed from the middle positions are shown with light-gray 
color. 



II. 1-DOPED-ELECTRON, 4-PLAQUETTE 
LINEAR CLUSTERS: BOND-LENGTH 
FLUCTUATIONS AND INTER-SITE CHARGE 
TRANSFER 

Multiconfiguration, complete active space (CAS) self- 
consistent-field (SCF) electronic structure calculations 
were carried out on linear clusters including four CUO4 
plaquettes. The CAS SCF method is designed to pro- 
vide an unbiased description of systems where near- 
degeneracy occurs among several electron configurations 
|20j . It is based on a partitioning of the total orbital space 
into three subsets, inactive, active, and virtual orbitals. 
The inactive orbitals are doubly occupied in all configu- 
rations and should be chosen as the orbitals that are not 
expected to contribute to near-degeneracy correlation ef- 
fects. The active orbitals are subject to no restriction 
on their occupation numbers and the virtual orbitals are 
unoccupied in all configurations. The CAS SCF wave- 
function is constructed as a full configuration-interaction 
(CI) expansion in the configurational space spanned by 
the active orbitals. Since both the CI coefficients and the 
orbitals are variationally optimized, the CAS SCF wave- 
function is extremely flexible. It is this flexibility that 
makes the CAS SCF method well suited for the descrip- 
tion of competing valence structures, bond breakings, 
and dissociation processes. 
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We used the crystal parameters reported for the 
electron-doped copper oxide superconductor compound 
Sri_ y NdyCu0 2 [23, with y w 0.16, Cu-0 distances of 
1.97 A, and no apical oxygens. The 4-plaquette clus- 
ter was embedded in an array of point charges that 
reproduce the Madelung potential associated with an 
idealized SrCuC>2 ionic system. We imposed Ci v sym- 
metry restrictions, with xy and yz mirror planes (see 
Fig. 2), and employed the following gaussian-type, atomic 
natural orbital basis sets for the [CU4O13] cluster: Cu 
21sl5pl0d6//5s4p3cfl/ and O 14s9p4d/4s3pld 0. It 
is known by now that such basis sets provide a quite ac- 
curate description of both ground-state and excited-state 
properties in a large variety of transition metal (TM) ox- 
ides, see for example H H |3 and 0. The Cu 
and Sr nei ghb oring ions were represented by total ion 
potentials |2(|. The calculations were performed with 
the MOLCAS program package |27| . 

With one electron added to the undoped, formally 
Cud 9 Op 6 configuration, the minimal active orbital space 
should consist of four (or even three) open-shell 3d x 2_ y 2 
orbitals. However, trial CAS SCF calculations were car- 
ried out with few virtual orbitals added to the active 
space. It turned out that these extra orbitals have either 
Cu 4s or Cu 3d' character, where 3d' is a more diffuse 
orbital than 3d [2^ , and non- negligible occupation num- 
bers. Also, large effects were observed on the absolute 
energies, as compared to the minimal active space. 

The relevance of the 4s orbital to the local-density 
approximation (LDA) conduction bands and to the in- 
plane/inter-plane hoppings was analyzed for several cop- 
per oxide com pou nds in Refs. HIES. A clear correlation 
was observed [3(| between the value of the critical tem- 
perature T c and the energy and specific composition of 
the so-called axial orbital, a hybrid with mixed Cu 4s, 
Cu 3e? 32 2_ r 2, and apical-oxygen 2p z character. The im- 
portance of the Cu 4s and Ap orbitals for an accurate 
description of the electron occupation of the Cu 3d x i_ y i 
and O 2p x /2p y orbitals in model Hamiltonian calcula- 
tions and also their role in screening the on-site Cu 3d 
Coulomb repulsion was discussed in 

We investigated the dependence of the cluster total 
energy on various distortions and varied the positions of 
the Cu ions and of those oxygens (Of,, see Fig. 2) bridg- 
ing the cations. For such a cluster, energetic effects as- 
sociated with displacements of the O a oxygen ions, see 
Fig. 2, are not meaningful, because part of the nearest 
neighbors of these ligands are modeled by effective po- 
tentials or even point charges. Nevertheless, under the 
assumption of charge segregation into filamentary stripe 
segments and one-dimensional conduction 0, 0, fl3| . 
the distortions along the y axis should be the most rele- 
vant. We found that the cluster energy is minimized for 
a geometry where the Cu-Of, bonds on one of the CUO4 
plaquettes are stretched by approximately 8% and the 
nearest Cu neighbors are also displaced by about 1% of 
the initial Cu-0 bond length, each of them in the same 
direction as the adjacent oxygen. The doped electron is 



accommodated in the 3d x 2_ y 2 orbital on the plaquette 
with elongated Cu-Ob bonds to give a Cu 3d 2 2 _ y2 con- 
figuration. We report results for geometries where cither 
the Cuz-Ot, or Cu r -0& bonds are stretched, see Fig. 2, and 
discuss the changes in the charge distribution for these 
two plaquettes and the variation of the cluster energy. 
However, in this study we keep all Cu ions at fixed po- 
sitions. This has only minor effects on the ground-state 
wave-function and energy. Our active space includes four 
3d x a_ y a, the Cu; 3d 3z 2_ r 2 and Cu r 3d 3z 2_ r 2 orbitals, and 
four virtuals. This is a 9 electrons/ 10 orbitals CAS, with 
about fifty thousand determinants. 

For the minimum-energy state with stretched Cu;-Of, 
bonds and a Cu; 3d 2 2y2 electron configuration, for ex- 
ample, those four correlating orbitals taken from the 
virtual orbital space in order to enlarge the minimal 
CAS have Cu; 3d! x2 _ y2l Cu; 4s, Cu r 3d' 3z2 _ r2 and Cu r 
3d' x2 _ y2 character. There is some mixing between the Cu; 
3c?3 Z 2_ r 2 and 4s basis functions. The occupation of the 
Cu; 3d 3z 2_ r 2 -like orbital is 1.97 and the predominantly 
Cu; 4s orbital has an occupation number of 0.03. There 
are also large contributions of the Cu 4s functions to 
some of the other active orbitals. The Cu r 3(^2 _ r 2 and 
3d' 3z 2_ r 2 orbitals have occupation numbers of 1.99 and 
0.01, respectively. The occupation numbers of the Cu; 
3d' x2 _ 2 and Cu r 3d' x2 _ y2 orbitals are somewhat smaller 
than 0.01. 

Inclusion of the diffuse d' orbitals into the active space 
is necessary for describing changes in the radial charge 
distribution for different d n configurations in the multi- 
configurational wave- function. For free TM atoms (ions) , 
this is a correlation effect involving pairs of electrons in 
the 4s 2 3d™, 4s3d" +1 , and 3d n+2 states @. Occupa- 
tion numbers as large as 0.1 for the d' orbitals of the Ni 
atom 31] demonstrate that a careful analysis is needed 
when dealing with TM atoms or TM compounds. In 
the copper oxide systems, the different d n configurations 
may be related to both on-site Cu 3d-4s and inter-site, 
Cu-Cu and O-Cu, excitations. Radial, or "breathing", 
correlation effects may also occur for the O 2p orbitals. 
When th e lig and to metal charge transfer excitations are 
relevant [32] , a highly accurate CAS SCF wave- function 
would require the inclusion of two extra sets of orbitals 
into the active space, ligand 2p orbitals involved in the 
2p-Cu3d bonding plus the corresponding 2p' orbitals 
[25II32]] . Such an extension of the active space is not pos- 
sible for our cluster, due to the high computational cost. 
Nevertheless, the low occupation numbers of the 3d' 2 , 2 

y 

orbitals, see the preceding paragraph, suggest that the 
essential near-degeneracy correlation effects are already 
included in our model. It will be shown in the next sec- 
tion that the role of the O 2p' orbitals becomes important 
in the hole-doped case. 

The potential energy landscape when stretching the 
Cu;-Ob or Cu r -Of, bonds is shown in Fig. 3 . The middle 
point on the horizontal axis corresponds to the undis- 
torted cluster. From each minimum-energy configura- 
tion, the two O;, ions on the Cu 3d 2 2 _ 2 plaquette are 
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FIG. 3: Results of static, total energy calculations for an 
electron-doped, 4-plaquette cluster, see text and Fig. 2. The 
filled symbols are the CAS SCF data. The lowest points 
correspond to Cuj 3d\_„2 or Cu r Sd\_„2 configurations 
on plaquettes with Cu-Of, distances larger by 8%. For ge- 
ometries where the Cu-Ot, bonds on the adjacent plaquette 
are stretched by more than 12%, an 1-electron Cu; d x 2_ y 2- 
Cu r d x 2_ y 2 charge transfer occurs. The dotted curves are 
obtained by CAS SI calculations. The CAS SCF energies for 
the structures with CuOt, bonds stretched by 8% on one of 
the Cu;04 or Cu r 04 plaquettes are taken as reference, a is 
the in-plane lattice constant. 



shifted back to the undistorted structure in steps of 2% 
of the high-symmetry Cu-0 distance. The orbitals from 
the previous step are used each time as starting or- 
bitals. Once the undistorted geometry is reached we start 
stretching the Cu-Of, bonds on the adjacent plaquette, 
see Fig. 2. Each curve in Fig. 3 is obtained by following 
this procedure. It turns out that the initial electronic 
wave-function remains largely unchanged for a relatively 
broad range of Cu-Of, distortions. Its character changes 
abruptly when the Cu-Of, bonds on the adjacent plaque- 
tte are elongated by about 12%. Only for such a geome- 
try the doped electron is transferred to the other d x 2_ y 2 
orbital. 

We list in Table I Mulliken populations of atomic or- 
bitals on the CU/O4 and Cu r 04 plaquettes for the Cuf 
d x \_ y 2 lowest energy configuration. The charges associ- 
ated with the O 2s and O 2p. x /2p y functions overlapping 
the Cu d x 2_ y 2 orbitals are sensibly smaller than 2.00, 
even on the d^ 2 _ y 2 plaquette. The "missing" charge is in 
the Cu 4s and 4p orbitals. This deviation from the fully 
ionic picture was previously pointed out in Ref. 0, for 
example. Between the two plaquettes, the major differ- 
ences are related to the populations of the Cu 3d x 2_ y 2, 
3rf 3z 2_ r 2, 4s, 4p, and O a 2p x orbitals. 

If the matrix elements between the Cu; d\_, 2 and Cu,. 

x y 

d 2 2 2 electronic wave- functions are ignored, the adia- 

x — y ° ' 

batic potential energy curves for the O-ion half-breathing 



displacements discussed above intersect at the point as- 
sociated to the undistorted cluster, see Fig. 3 . However, 
if the electronic structures are allowed to relax by mix- 
ing of the Cu; d^2_ v 2 and Cu r d^ 2 _ v 2 configurations, 
no intersection would occur. A so-called avoided cross- 
ing is expected for the potential curves. The matrix 
elements between non-orthogonal, separately optimized 
CAS SCF wave- functions can be obtained by State In- 
teraction (SI) calculations [13 ■ The potential energy 
curves corresponding to the CAS SI states are shown in 
Fig. 3 by dashed lines. The interaction between the non- 
orthogonal states is not too large, with energetic effects of 
not more than few tenths of eV. The weights of the non- 
orthogonal states in the lowest CAS SI root are both 1/2 
at the middle point. For chain bonds stretched by 8% on 
one of the plaquettes, the weights are 0.99 and 0.01 . The 
energy barrier associated with the lower CAS SI curve is 
540 meV. 

For a more accurate estimate of the size of the energy 
barrier, polarization effects beyond the boundaries of the 
linear, 4-plaquette cluster should also be taken into ac- 
count. Test calculations on a 5-plaquette cluster where 
all four nearest neighbor plaquettes of a "central" Cu c 
d^2_,2 ion are treated at the all-electron level show in- 

x y 

deed that such effects are important. We used for these 
calculations D%h symmetry restrictions and the following 
basis sets 22]: Cu 21sl5pl0d6//5s4p3d2/ for the five 
TM ions, O 14s9p4d/4s3p2d for the ligands bridging two 
cations and O 14s9p/4s3p for the other twelve oxygens. 
For a Cu c d^ 2 _ y 2 configuration on the plaquette at the 
center of the cluster, the CAS SCF energy difference be- 
tween the undistorted structure and the geometry where 
the Cu c -0 distances along the y axis are longer by 8% 
is 130 meV lower than for the linear, 4-plaquette cluster. 
Applying such a correction to the CAS SI results, the en- 
ergy barrier would decrease to about 400 meV. Longer 
range polarization effects should further lower this value. 

We also performed calculations on a 6-plaquette lin- 
ear cluster with two doped electrons. The following 
contractions were applied for the Cu and O basis sets: 
Cu 21sl5pl0d/5s4p3d and O 14s9p/4s3p 22]. Ground- 
state Mulliken populations for the geometry with Cu-Of, 
bonds stretched by 8% on the Cuzz04 and Cu rr 04 pla- 
quettes are given in Table II. The results were obtained 
with an active space including 12 electrons and 12 or- 
bitals. For this geometry and choice of the multiconfigu- 
rational space, the CAS SCF active orbitals were identi- 
fied as six Cu 3d x 2_ y 2 orbitals, the Cuu 3c?3 Z 2_ r 2, 4s and 
Cu rr 3c?3 Z 2_ r 2, 4s plus two correlating Cuu 3d' x2 2 and 
Cu rr 3g?'j, 2 . The occupation numbers of these active 
orbitals are similar to those found for the smaller cluster. 
The Mulliken populations are very similar, too. The re- 
sults show thus that localized 3d%_ 2 states are also ob- 

x y 

tained in larger clusters and for more than a single doped 
electron. For the electron-doped case, an analysis of 
inter-carrier couplings by CAS SCF calculations on such 
larger clusters will be provided in a forthcoming publi- 
cation. In the case of hole-doping, preliminary CAS SCF 
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TABLE I: Mulliken charges per centre and basis function type for the minimum-energy Cu; d^ 2 _ y 2 configuration on an electron- 
doped, 4-plaquette linear cluster. The Cuj-Of, bonds are stretched by 8%. Values for the Cu;04 and Cu r 04 plaquettes are 
listed. The Mulliken populations on the other plaquettes, Cu«04 and Cu r r04, are very similar to those on Cu,-04. The charge 
of each 3d t2 g atomic orbital is 2.0 . There is also some charge associated with the Cu / and O d functions, not shown in the 
table. 
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TABLE II: Mulliken charges per centre and basis function type for a 6-plaquette linear cluster with two doped electrons. The 
Cujj-Of, and Cu rr -0(, bonds are stretched by 8%. Values for Cu and Of, ions on the four central plaquettes, Cu„04, Cu,04, 
Cu r 04, and Cu rr 04, are listed. The same notations as in Fig. 2. are used. 
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calculations on 2-oxygen-hole, 6-octahedra clusters re- 
veal that the hole-hole interaction is large. Due to the 
strong mutual repulsion within the O 2p "band" , the 
two holes are pushed to the boundaries of the cluster 
and meaningful CAS SCF results cannot be obtained in 
this situation. However, inter-carrier correlations for the 
hole-doped case are addressed in the next section by a 
model Hamiltonian approach that enables calculations 
on chains of few tens of plaquettes. 

The existence of a local double-well potential like 
that shown in Fig. 3 for Cu-0 bond-length distortions 
in both n-doped and p-doped layered Cu oxides was pre- 
dicted by Goodenough [l3|- Goodenough also suggested 
that around optimal doping the doped carriers display 
some particular (dynamic) ordering within each layer, 
perhaps an arrangement with alternating undoped and 
1/3-doped chains of plaquettes. The coupling between 
these ordered electrons/holes and in-plane O-ion half- 
breathing phonons is believed to be at the heart of high- 
temperature superconductivity [Til Il2l Il3| . Evidence for 
stripe formation in electron-doped cuprates was recently 
reported in Ref . |34j . for example. 



III. HOLE-DOPED SYSTEM 

Multiconfiguration calculations for hole-doped clusters 
have been performed in Ref. 6j. However, a minimal ac- 
tive space was applied in that study, with a single orbital 



for each Cu 3d or O 2p hole. In this section we report re- 
sults of CAS SCF calculations with larger active spaces. 
We use the same 4-octahedra linear cluster and the same 
(static) embedding as in R ef. B , with the lattice param- 
eters of Lai.85Sri. 85 Cu04 [33. The basis sets are also 
the same, Cu 21sl5pl0d/5s4p3d, O Us9p4d / 4s3pld for 
the bridging Ob oxygens, and O I4s9p/4s3p for the other 
oxygens |22j. As already mentioned in Introduction, the 
CAS SCF results indicate that ZR type singlet states 
are only formed on CUO4 plaquettes with Cu-0 bonds 
shorter than the average Cu-0 distance deduced from 
experiment. On the other hand, for the undistorted ge- 
ometry the doped hole has the largest weight onto a sin- 
gle ligand and in the simplest picture the dominant con- 
tribution to the wave- function comes from a ...-Cud 9 - 
Op 5 -Cud 9 -... configuration (3(|. The latter state is few 
hundreds meV higher in energy than the former. Also 
for the hole-doped case, O-ion half-breathing displace- 
ments determine thus a multi-well energy landscape and 
can induce charge transfer along a chain of CUO4 plaque- 
ttes. We note that geometry optimizations for transition 
metal oxides like MnO and NiO at the Hartree-Fock (HF) 
level [13, HU slightly overestimate the interatomic dis- 
tances. This means that the multi-well energy landscape 
would also be obtained with such HF optimized lattice 
parameters 39]. For the electron-doped case, the usual 
overestimation of the lattice constants, 2 — 4% |3H l38j. 
is significantly lower than the range of the Cu-0 bond 
distortions associated with the multi-well energy curve, 
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TABLE III: Major contributions of the oxygen 2p and metal 3d atomic orbitals to those active (natural) orbitals implying 
2p—2p and 2p— 3d bonding on the Cu;04 and Cu r 04 plaquettes. CASSCF results for a hole-doped, 4-plaquette linear cluster. 
Coefficients smaller than 0.20 are not listed. For each geometry, the nature of the other active orbitals is discussed in the text. 
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Oj-hole, undistorted cluster 






Cu;04 ZR, distorted Cuf-Ob bonds a 




1.99 


1.96 


1.84 


0.21 


1.99 


1.97 


1.90 


0.13 


O l b Py 




0.54 


0.22 


0.37 


0.29 


0.67 


0.41 


0.42 


O l a p x (x2) 


±0.20 


±0.25 


±0.26 


±0.31 


T0.42 




±0.29 


±0.29 


Cuj d x 2 y2 






-0.61 


0.74 






-0.63 


0.77 


Ot Py 


0.37 


0.51 


-0.56 


-0.47 


-0.29 


0.67 


-0.41 


-0.42 


Cllr d x 2 y 2 








-0.21 










o; p, (x2) 


=F0.37 


±0.21 














0£ Py 


0.39 
















a The chain Cu ; 


-Of, bonds arc 


6% shorter. 















±8%. 

As in the previous section, we added to the minimal 
active space orbitals from both the inactive, doubly occu- 
pied set and the virtual space. Trial calculations showed 
that the most important near-degeneracy correlation ef- 
fects are already accounted for with extensions of only 
few orbitals. In Tables III, IV, and V we show results ob- 
tained from CAS SCF calculations with 13 electrons and 
12 orbitals in the active space. Compared to the minimal 
(i.e., 5 electrons/5 orbitals) CAS, four inactive and three 
virtual orbitals were moved into the active set. Depend- 
ing on the input starting orbitals, different CAS SCF so- 
lutions can be obtained for the undistorted cluster, with 
the doped hole having the largest weight onto one of the 
three bridging oxygens, O l b , or 0£. The largest con- 
tributions to those active orbitals implying 2p-2p and 
3d-2p bonding on the central Cuf04 and Cu r 04 pla- 
quettes for the CAS SCF solution with the doped hole 
having the largest weight onto the ligand are illus- 
trated in Table III. We also display in Table III the oc- 
cupation numbers of those orbitals. There are other 
three active orbitals with Cuzf, Cu r , and Cu rr 3d x 2_ y 2 
character and occupation numbers of nearly 1. For this 
CAS SCF solution, the remaining active orbitals have O b 
2s, Cu; 3d 3z 2_ r 2, O c b 2p' y: Cuf4s-0£2s'-Cu r 4s, and 
Cu; 3d' 3z 2_ r 2 character. Their occupation numbers are 
1.999, 1.98, 0.02, 0.01, and 0.01, respectively. 

The Mullikcn electron populations for the relevant 
atomic functions on the Cu;04 and Cu r 04 plaquettes are 
listed in Table IV. Compared to the minimal-CAS calcu- 
lations from Ref. , the charge disproportion between 
the (O l b 2p y , O l a 2p x ) and 2p y orbitals is significantly 
smaller. The difference between the O l b 2p y and O b 2p y 
Mulliken charges, for example, is 0.13 instead of 0.50. 
We also note the Cu;- O b - Cu r asymmetry. A symmetry 
equivalent solution, with a pair of bonding and antibond- 
ing O b 2p y - Cu r 3d x 2_ y 2-... orbitals that is symmetry 
equivalent to the pair of bonding and antibonding 
Cu; 3d x 2_ y 2- O b 2p y -... orbitals in Tablelll (the third 
and fourth columns) can also be obtained. The CAS SCF 
wave-function changes its character from one type to the 
other for displacements of the 0£ ion of ±0.8% around the 



middle position. This is very different in comparison with 
the electron-doped case in the previous section, where the 
character of the CAS SCF wave-function remains nearly 
unchanged, i.e. either Cu; d^ 2 _ y2 or Cu r d^ 2 _ y2 , for dis- 
tortions of up to ±12% of the Cu-Oft bonds. 

With the larger active space, the ZR like state with 
shorter Cu-0& distances (6% shorter) on one of the 
Cu;04 or Cu r 04 plaquettes is lower in energy by about 
460 meV. Major contributions to the relevant 2p—3d ac- 
tive orbitals are shown in Table III. The doped hole is now 
approximately equally distributed over the Of, p y and O a 
p x atomic orbitals on the distorted plaquette, see the pair 
of d—p bonding and antibonding orbitals with lower occu- 
pation numbers in Table III. As above, there are three ac- 
tive orbitals with Cu;;, Cu r , and Cu rr ?>d x 2_ y 2 character 
and occupation numbers of roughly 1. The other active 
orbitals have Cu; 3p y , Cu/ 3d Zz2 _ r2 , Cu; 4s — O b / a 2s'/2p', 
O^, 2p' y -Cui 4p y -Og 2p y , and Cu; 3d' 3z2 _ r2 character and 
occupation numbers of 1.999, 1.98, 0.01, 0.01, and 0.01, 
respectively. It can be seen that the nature of some of 
these orbitals changes from the undistorted to the dis- 
torted cluster. The CASSCF results indicate that in 
the latter case the 2s/2p-Cu4s/4p-0 2s/2p interac- 
tions on the same plaquette are more important than 
the inter-plaquette Cu4s/4p-0 2s/2p-Cu4s/4p inter- 
actions for the former. Mulliken charges associated with 
such a ZR like state are listed in Table V. The compari- 
son between the results displayed in Table V (or Table IV) 
and those obtained for the electron-doped CU/O4 plaque- 
tte in Table I reveals a remarkable feature: due to strong 
charge redistribution effects within the Cu/ 3d 3z 2_ r 2, 4s, 
4p and O 2s, 2p orbitals, the total electronic charge asso- 
ciated with the Cu; Sd^. 2 _ y2 ion in Table V (or Table IV) 
is nearly as high as that associated with the Cu; 3d^ 2 _ y 2 
configuration from Table I, 28.7 versus 28.8. It is well 
known that the Mulliken charges may be misleading in 
cases where strong inter-atomic orbital overlap occurs. 
Still, a variation of nearly one elementary charge within 
the Cu 4s, 4p shell cannot be an artifact. 

The adiabatic energy landscape for Cu-O;, bond-length 
distortions on the two central plaquettes is shown in 
Fig. 4. The points on the full lines are the minimal-CAS 
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TABLE IV: Mulliken charges per centre and basis function type for a broken-symmetry oxygen-hole state on a 4-octahedra 
linear cluster. The doped hole has the largest weight on the 0£ ion. Charges for the central CuiCU and Cu r 04 plaquettes are 
listed. The Cu-O distances are all the same. 





oi 


oi 


Cu; 


ot 


Cu r 




01 


O 2s, Cu 4s 


1.59 


1.81 


0.77 


1.58 


0.67 


1.81 


1.60 


O 2p x , Cu 4p x 


2.00 


1.65 


0.28 


2.00 


0.26 


1.81 


2.00 


2p y , Cu 4p y 


1.65 


2.00 


0.22 


1.51 


0.17 


2.00 


1.83 


2p z , Cu 4p z 


1.93 


1.97 


0.25 


1.92 


0.23 


1.97 


1.93 


Cll d$ z 2 _ r 2 






2.00 




2.00 






Cu d x 2_ y 2 






1.06 




1.15 






Total el. charge 


9.19 


9.42 


28.70 


9.03 


28.57 


9.59 


9.37 



TABLE V: Mulliken charges per centre and basis function for a minimum-energy ZR type configuration on the Cu;04 plaquette. 
CAS SCF results for an 1-doped-hole, 4-octahedra linear cluster, see text. The Cui-Oi, bonds are shorter by 6%. Values for 
the central Cu(04 and Cu r 04 plaquettes are listed. 





o l b 


O l 


Cu ; 


Ol 


Cu r 




o r b 


O 2s, Cu 4s 


1.57 


1.81 


0.82 


1.57 


0.64 


1.80 


1.60 


O 2p x , Cu 4p x 


2.00 


1.65 


0.29 


2.00 


0.25 


1.83 


2.00 


O 2p y , Cu 4p y 


1.59 


2.00 


0.25 


1.59 


0.15 


2.00 


1.84 


O 2p z , Cu 4p z 


1.92 


1.97 


0.25 


1.92 


0.23 


1.97 


1.93 


CU d$ z 2 _ r 2 






2.00 




2.00 






Cu d x 2_ y 2 






1.02 




1.16 






Total el. charge 


9.10 


9.43 


28.74 


9.09 


28.53 


9.61 


9.38 



results already published in Ref. [fj. The energy max- 
ima correspond to states with dominant Ob-hole charac- 
ter in the undistorted structure. The two minima are 
related to ZR like states where the Cu-0& distances on 
the CU/O4 plaquette (the left hand minimum) or the Cu- 
Ob distances on the Cu r 04 plaquette (the right hand 
minimum) are 6% shorter. From each minimum-energy 
structural configuration, the two Of, ligands are shifted 
back to the undistorted geometry in steps of 1% of the 
high-symmetry Cu-0 bond length. These states are 
shown with black dots. The last point on each of the 
two curves corresponds to a configuration where the Cu- 
Ob distances on either the Cu;;04 or Cu rr 04 plaquette 
are shorter, by 1%. Other calculations were carried out 
for the case where the Cu-Ob bonds on the CU/O4 and 
Cu r 04 plaquettes are gradually shortened starting from 
the undistorted structure and using the orbitals of the 
Oj-hole state as starting orbitals. Due to convergence 
problems, only few points could be obtained. Those 
points are shown as black triangles. 

Relative energies from CAS SCF calculations with 13 
electrons and 12 orbitals in the active space are shown 
in Fig. 4 only for the ZR like states with Cu-O;, distances 
shorter by 6% and the states with undistorted bonds. 
These points are marked with filled squares. The energy 
lowering induced by the use of a larger active space is of 
the order of 1.5-2.0eV. The effect is stronger for the dis- 
torted configurations. For convenience, the larger-CAS 
energies are shifted upwards in Fig. 4 by 1.6 eV. We note 
that the active orbitals and the CI vectors for the states 
with dominant O^ or 0£ -hole character are similar to 



those of the state where the hole has the largest weight 
on the Og ligand, with a "shift" of one plaquette to the 
left or to the right, respectively. 

Separately optimized, broken-symmetry wave- 
functions like the different O-hole solutions for the 
undistorted geometry are non-orthogonal and interact- 
ing. Orthogonal, non-interacting eigenstates can be 
obtained by SI calculations. However, the situation is 
now somewhat different as compared to the SI calcu- 
lations for the electron-doped cluster. First, the way 
the energy landscape was drawn is different. The two 
sides of a well in the hole-doped case correspond to 
identical positions of the Ob ions, although CAS SCF 
roots on different sides have different character, at least 
for small distortions. In the second place, the O-hole 
is "smoothly" transferred from one plaquette to the 
other, for distortions of about 1% of the Cu-O bond(s). 
Matrix elements exist thus only between different O-hole 
states in the undistorted geometry or between states 
situated on different sides of the same well. We note 
that in the electron-doped cluster, matrix elements 
exist only between CAS SCF states associated with 
different wells. Those matrix elements determine in 
the electron-doped case a so-called avoided crossing 
for the energy curves. The missing points for each 
of the two wells in Fig. 4 could probably be obtained 
from calculations on clusters including several 1-hole, 
3-octahedra units. For certain phases of the O-ion 
half-breathing vibrations on the neighboring plaquettes, 
see Fig.l, the O-hole on a given 3-octahedra unit should 
move "smoothly" from one side to the other. In fact, 
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it is the longer range Coulomb and spin interactions 
that can induce indeed broken-symmetry oxygen-hole 
states in the infinite system. The collective, coherent 
"tunneling" of these quasi-localized oxygen holes along 
a multi-well energy landscape like that depicted in Fig. 4 
can lead to a traveling charge-density/spin-density wave 
as in the model proposed by Goodenough [l2| . 

The tendency towards the formation of quasi-localized 
oxygen-hole states on distorted CUO4 plaquettes is con- 
firmed by model Hamiltonian calculations. We consider 
in the following paragraphs isolated chains of CUO4 pla- 
quettes at 1/3 hole-doping and construct &p—d, Hubbard 
type model: 

H = e d ^2 n i<y + e v X! n i° + X] ( d i<yP^ + h ' c -) 

<jj'><T 

i 3 

T Vpd ^ ^ 'flifjTljfj' -)- Vpp ^ ^ TljfjTljifjt . 
<ij>aa f <jj'>aa' 

We neglect the apical ligands and restrict to a single 
orbital at each site: Cu 3d x 2_ y 2 , bridging oxygen Of, 
2p y orbitals, and oxygen 2p x orbitals at the other lig- 



500 




-200 - ■ ■ 

-6(1) -6(r) 

[2xd(CU| r O b )-a]/a(%) 



FIG. 4: Results of static, total energy calculations for a 
hole-doped, 4-plaquette cluster, see text. Filled circles and 
triangles are the minimal-CAS (i.e. 5 electrons/5 orbitals) re- 
sults. The lowest points are related to ZR like states where 
the Cuj-Oi, distances (the left hand minimum) or the Cu r - 
Ot distances (the right hand minimum) are 6% shorter. The 
energy maxima correspond to broken-symmetry states in the 
undistorted geometry. The filled squares are CAS SCF re- 
sults with 13 electrons and 12 orbitals in the active space. 
The minimal-CAS energy of the ZR state is taken as refer- 
ence. For convenience, the larger-CAS energies are shifted 
upwards by 1.6 eV. 



and sites, O a . We take into account the overlap integrals 
and Coulomb interactions between the nearest neighbor 
d x 2_ y 2 and p y /p x er-orbitals and also between the sec- 
ond nearest neighbors Ob p y and O a p x . In the electron 
representation, we choose A = e d — e p = 3.0 [29| and set 
t x pd = -1.2, t y pd = -1.2, and t pp = 0.5 (eV). The values 
of the Coulomb repulsion parameters are chosen accord- 
ing to results of constrained LDA |40| and/or ab initio, 
explicitly correlated 0, 0, 0, calculations: U d d = 9.0, 
Upp = 5.0, V p d = 1-5, and V pp = 1.2 (eV). Under the as- 
sumption of c harg e segregation into undoped and 1/3- 
doped stripes [12], we a l so include for the Q ligands 
inter-site Coulomb interactions due to electrons at Cu 
and O sites on the neighboring, undoped Cud 9 - Op 6 -... 
chains. That is, the energy of the O a p x orbitals is shifted 
from e p = e^ = e d -3 to = e p + V pd + 4 V pp . 

Ground-state expectation values are obtained by 
density-matrix renormalization-group (DMRG) tech- 
niques 41]. For these calculations chains of 33 CUO4 
plaquettes, open-end boundary conditions, and to =1200 
states to build the DMRG basis were used. The typical 
discarded weight is of the order 10~ n — 10~ 10 or lower. 

Occupation numbers of p and d orbitals on the 3- 
plaquette unit at the middle of the 33-plaquette chain 
are listed in Table VI (first column). The Cu site right at 
the middle is labeled Cu c . Along the Cu-O;,-... chain, 
the next sites to the left and to the right are denoted 
O l b , Cu;, O" and 0£, Cu r , 0^ r , respectively. For each 
of the Cu ions mentioned above there is also a pair of 
O a sites. These ligands are labeled O^, O^ or O^. The 
occupation numbers of the d orbitals are comparable to 
those obtained in the CAS SCF calculations for the small 
cluster. However, since the O 2s, 2p — > Cu4s, Ap (back) 
charge transfer mentioned above is not taken into account 
in our p— d model Hamiltonian, the electron occupation of 
the p orbitals is somewhat larger. Also, compared to the 
ab initio results, there is no charge disproportion along 
the undistorted chain of CUO4 plaquettes. This may be 
attributed to various effects. The most important should 
be the neglect of longer range Coulomb interactions in 
the p — d model Hamiltonian and the orbital optimiza- 
tion for the CAS SCF wave-function. The fact that the 
variational optimization of the "molecular" orbitals can 
lead to broken-symmetry solutions was shown for a sim- 
ple model Hamiltonian by Martin |4|. 

Next, we simulate the shortening of the Cu-O;, bonds 
on every third plaquette along the chain by introducing a 
modulation of the hopping matrix elements and inter-site 
Coulomb repulsion parameters: t' pd (r) =t p d(0) (r /r) 7 / 2 , 
t'pp{r)=t P p{0){r o /rY 0, and V (r) = V(0) (r„/r). Or- 
bital occupation numbers are given in Table VI, for bonds 
shortened by 6% on the central plaquette on every 3- 
plaquette unit (second column) and for bonds shortened 
by 6% on the right-hand plaquette on every 3-plaquette 
unit (third column). We also list the energies calculated 
for the three configurations, with the energy of the undis- 
torted chain taken as reference. The results show that the 
energy of the system is lowered for the distorted config- 
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TABLE VI: Orbital occupation numbers on the 3-plaquette 
unit at the middle of a chain of 33 CUO4 plaquettes at 1/3 
hole doping. DMRG results for a p — d extended Hubbard 
model, see text. Three different geometries are considered. 
In the distorted configurations, the Cu-Oj, bonds are shorter 
by 6% on every third plaquette. Each occupation number 
does not change by more than 0.02 for the next four "units". 
Relative energies are also given, for a whole chain. 





Undist. 


Dist. Cu c -Oi, 


Dist. Cu r -0(, 


O'^Py 


1.75 


1.80 


1.71 


Cu; d x 2_ y 2 


1.14 


1.18 


1.18 


oi Px 


1.89 


1.93 


1.93 


0{ Py 


1.74 


1.71 


1.81 


Cu c d x 2_ y 2 


1.13 


1.04 


1.18 


01 Px 


1.88 


1.83 


1.93 


01 P y 


1.74 


1.71 


1.71 


CUr d x 2 y 2 


1.14 


1.18 


1.04 


O r a P X 


1.89 


1.93 


1.83 


or P y 


1.75 


1.80 


1.71 


Rel. En. (eV) 





-2.80 


-2.44 



urations. The small energy difference between the two 
differently distorted structures, 2.80 versus 2.44 eV, is 
due to finite size effects. In addition, the charge redistri- 
bution within the 3-plaquette unit(s) confirms the trend 
observed for the CAS SCF results on the smaller cluster: 
the oxygen holes tend to localize on the distorted pla- 
quettes. Electron-phonon interactions and a softening of 
the half-breathing O-ion vibrations were predicted in the 
hole-doped regime by other authors as well ^3|> mostly 
on the basis of t— J model Hamiltonian studies. With the 
present DMRG calculations for a p—d extended Hubbard 
model, we found that the energetic effects associated with 
Cu-Oh bond-length distortions on every third plaquette 
along a chain of several plaquettes at 1/3-doping is of 
the order of hundred meV per each 3-plaquette "unit" 
(approximately 255 meV for distortions of 6%). These 
estimates are in line with the CAS SCF results discussed 
above. Nevertheless, the ab initio calculations indicate a 
much stronger effect. 

IV. CONCLUSIONS 

Results of ab initio multiconfiguration calculations on 
clusters including few CUO4 plaquettes and also results of 
DMRG calculations for a p—d extended Hubbard model 
with chains of few tens of plaquettes indicate the presence 
of strong anharmonic effects for the O-ion half-breathing 
deformations. In addition, for certain phases of the O- 
ion vibrations, charge can be displaced along the CuO 
chain. These effects are observed for both hole and elec- 
tron doping. 

The conditions under which the Cu02 layers can 
be doped were discussed by Goodenough et ai, see 



Refs. 0, for example. Electrons are not easily ac- 
cepted if the in-plane lattice parameter is too small, 
a < 3. 93 A. On the other hand, holes can be introduced if 
the lattice parameter is not larger than 3.87 A. Strong 
couplings between the Cu-0 bond-length fluctuations 
and the charge carriers, with an anharmonic, multi-well 
energy landscape and strong charge transfer effects were 
also predicted in Ref. |12J. This prediction is confirmed 
by our present study. Nevertheless, both the CAS SCF 
results on the small clusters and the DMRG estimates for 
the p—d, Hubbard like CUO4- CUO4-... chain are rather 
qualitative. For the CAS SCF calculations, a large part of 
the dynamic electron correlation effects is not accounted 
for jig. Coulomb and (super)exchange interactions be- 
yond the 4-plaquette cluster are also neglected. Within 
the hole-doped, p — d extended Hubbard model we have 
neglected "rehybridization" effects Q arising from the 
variational optimization of the orbitals. Charge dispro- 
portion along the CuO chain can be induced then only 
by longer range Coulomb interactions and structural dis- 
tortions. However, since interactions beyond the near- 
est neighbors are not explicitly included, longer range 
Coulomb "correlations" occur through indirect, higher 
order effects. Each of the two approaches has thus its 
limitations, but to some extent the two types of methods 
are complementary. Most importantly, they both predict 
that the electron - lattice couplings in the doped Cu02 
planes are strong and anomalous. The superconductiv- 
ity might be related to the collective, coherent tunneling 
of the doped carriers along the multi-well energy land- 
scape. Tunneling effects associated with an electron hole 
in a double-well potential were analyzed on an oversim- 
plified, hole-doped [Cu02]~ cluster model by Bishop et 
ai, see and references therein. The height of the en- 
ergy barrier was 400 meV in their model, similar to our 
results. The amplitude of the O-ion displacements was 
about two times smaller. 

Within the traveling CDW/SDW scenario illustrated 
schematically in Fig.l, the rigidity of the superconduct- 
ing condensate is ensured by the inter-carrier Coulomb 
repulsion and the relatively large AFM interactions be- 
tween the intervening Cu d 9 ions. For a single chain, 
these couplings are weakened for doping levels below 1/3. 
In the case of hole doping, it seems that at higher con- 
centrations an homogeneous metallic state is realized, like 
the 1/2-doped metallic chains in Lai.4sNdo.4oSro.i2Cu04 
p9| . The charge segregation into 1 /3-doped and undoped 
stripes at optimal dopin g 1121 i s supported by inelastic 
neutron scattering data [5(J,|5l|]. Nevertheless, the exis- 
tence of this type of charge ordering in the Cu02 plane 
is not definitely established. 

We thank A. Yamasaki, G. Stollhoff, G. Kalosakas, 
and A. T. Filip for fruitful discussions. L. H. acknowl- 
edges partial financial support from the Alexander von 
Humboldt Foundation. 
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